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Abstract 

We develop an approximation scheme for our worldsheet model of the sum of planar diagrams 
based on mean field theory. At finite coupling the mean field equations show a weak coupling 
solution that resembles the perturbative diagrams and a strong coupling solution that seems to 
represent a tensionless soup of field quanta. With a certain amount of fine-tuning, we find a 
solution of the mean field equations that seems to support string formation. 
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1 Introduction 



Last year we initiated a program (T|, for summing the planar graphs of large N c matrix field theories, 
based on a description of each planar diagram as a light-cone worldsheet path integral. Once this 
worldsheet representation was obtained, we could then realize the sum over planar graphs by 
coupling the worldsheet matter and ghost fields to a two dimensional Ising spin system. Summing 
over all spin configurations accomplishes the sum over all planar diagrams. One can regard the 
resulting two dimensional system as a noninteracting string moving in a background described by 
the Ising spin system. In this way our proposal arrives at a result suggested by the Maldacena 
conjecture ||, namely the association of a free string theory in a non-trivial background with the 
sum of planar diagrams. But instead of relying on open-string/closed-string duality to provide 
equations for the background, we are able to "read off" the background directly from the planar 
diagrams being summed. Thus we hope our approach represents an attack on the problem of large 
N c QCD complementary to the well studied AdS/CFT duality. 

The worldsheet construction of Q was restricted to 3 scalar field theory, but the methods 
were soon extended to the important case of pure Yang-Mills field theory in ||. The extension to 
supersymmetric field theories remains to be done. But here we return to the simple </> 3 theory as a 
useful arena for developing methods for extracting the physical properties of the worldsheet system, 
and hence those of the large N c limit. In particular we wish to apply the mean field approximation 
to the Ising spin system in order to get a qualitative understanding of the basic physics described by 
the sum of planar graphs. Such an approach has previously been applied in an attempt to extract 
nonperturbative information from lattice string theory Q in Ref. ||. 

At the outset we should take note of the fundamental instability of (p 3 theory following from the 
unboundedness of the cubic potential. As long as the free fields we perturb around are not tachyonic, 
this instability does not obstruct the evaluation of any Feynman graph, nor does it prevent one from 
contemplating graph summation Q. Although our real interest is summing the planar diagrams 
of a stable non-abelian gauge theory, we find the simplicity of 3 diagrams particularly valuable as 
a toy system for developing and testing our methodology. In particular, we can ask whether our 
methods are at all sensitive to the instabilities we know are there. There is even the possibility 
that a metastable vacuum in the finite N c case becomes stable in the N c — > oo limit If so, the 
planar sum might give a physically meaningful answer, and show instability only at higher order 
in the 1/N C expansion. 

We develop the mean field method in stages. After a brief review of the worldsheet formalism 
in Section 2, we begin in Section 3 with a treatment motivated by the formal continuum limit 
of the worldsheet model. The qualitative continuum physics is reviewed and the mean field is 
introduced. The basic idea is that the three worldsheet systems in the construction, namely the 
target space variables q the ghosts b, c and the Ising spins s, are first exactly solved in the presence 
of a homogeneous mean field. Then the value of the mean field is determined by minimizing the 
sums of the energies of the three systems. Because the mean field (f) is supposed to represent a 
fundamental variable that only takes on the values or 1, there is some freedom in the way it 
can enter the dynamics, when < 4> < 1, in the first approximation to the actions that describe 
the target space and ghost systems. Postponing this issue to Section 4, we instead use simple 
scaling arguments to fix a reasonable interpolation for the mean ghost +matter energy between its 
values for = and <f> = 1. The resulting mean field equation has a solution that supports string 
formation for a certain choice of parameters. For this choice there is a minimum of the energy that 
is associated with a finite non-zero string tension. 

In Section 4, we develop the mean field approach directly from the lattice system that was 
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the foundation of our mapping of planar diagrams to worldsheets. This treatment is thus more 
transparently related to the original field theoretic perturbation theory. Here we identify the mean 
field approximation as a certain saddle-point evaluation of the path integral. In this context, the 
ambiguities mentioned in the previous paragraphs are simply rearrangements of the path integration 
that are valid in the exact integral but lead to differences in the saddle-point approximation. Since 
the saddle-point approximation can be systematically corrected, there is an objective criterion for 
choosing the zeroth order action, namely the one for which the corrections are as small as possible. 
We base our choice in this section on how well it works at zero coupling, where exact calculations 
can be made. It turns out that the consequences of this choice are qualitatively very similar to the 
results of Section 3. The main difference can be traced to the (p 3 ^ 2 dependence of the energy in (|5l|), 
as opposed to the linear dependence on cj> in (|28l). This is due to the ambiguities in the application 
of the saddle point method mentioned above, and it does not lead to any qualitative change in the 
final results. Again there is a range of parameters, where the coupling g = 0(e) where the effective 
string tension is finite and non-zero. Here e is a temporary infrared cutoff that must be sent to 0. 
Although the coupling is tending to 0, this is not the perturbative regime, which is obtained by 
expanding about g = at fixed e and lattice cutoffs, and only removing the cutoffs order by order 
at the end of the calculation. For g/e — > the effective string tension becomes infinite, signaling 
the return to the perturbative region of weakly coupled field quanta. 

For g = 0(1) as e —* 0, the effective tension in our approximation is 0: the system is some 
kind of tensionless soup. We believe that the unphysical nature of this regime is a reflection of the 
inherent instability of the initial </> 3 theory. We find the existence of a regime with string formation 
interesting, even though it seems to require a certain amount of fine-tuning. 

Concluding remarks are given in Section 5. We also include two appendices. In the first we 
explain the tools for evaluating the Gaussian path integrals encountered throughout the article. In 
the second, we discuss a simple extension of the mean field description to slowly varying fields. 



2 Review of the Worldsheet Formalism for Field Theory 

The light-front components of any Minkowski vector will be written {v + ,v~ , v) or {v + , v~ , v k ). 
Here = (v° ± v 3 )/y/2, and the remaining components label the transverse directions. The 
Lorentz invariant scalar product of two four vectors v, w is written v ■ w = v • w — v + w~ — v~w + . 
We shall select x + to be our quantum evolution parameter, and we recall that the Hamiltonian 
conjugate to this time is p~ . A massless on-shell particle thus has the "energy" p~ = p 2 /2p + . 

The starting point for the worldsheet construction jl[ is the path integral representation of the 
(imaginary time) light-cone evolution operator for a free particle or field quantum 

J DcDbDq e- So (1) 

where the prime denotes d/da, and where Dirichlet boundary conditions q(0, r) = qo, q(p + ,7") = 
QMj b = c = are imposed on the world sheet fields. The target space coordinates q(tr, r) are 
related to the transverse momentum carried by the system by p = f£ daq' = qjv/ — qo, and 
the Dirichlet boundary conditions on q ensure the conservation of total momentum. The ghosts 
6, c, which are Grassmann variables, are necessary to ensure the correct measure factors. We shall 
always understand path integrals to be the continuum limit of ordinary integrals over variables 
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defined on a lattice We specify the lattice spacing in r to be a and that in a to be m. The 
continuum limit will be a, m —* with m/a fixed. Then the measure of the path integral is given 
by 
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where r = Na and p + = Mm, with M, N large positive integers. Note that since a has dimensions 
of momentum and r has dimensions of time, m/a has the dimensions of force. 

As discussed in a general planar diagram in quantum field theory can be represented as 
a path integral similar to (||) but with any number of internal Dirichlet boundaries given by any 
number of parallel line segments at fixed a, summed over varying location and length. Cubic 
vertices are simply represented by the appearance or disappearance of a Dirichlet boundary, and 
so are characterized locally on the world sheet. Vertices of higher order are represented by the 
simultaneous appearance and disappearance of more than one Dirichlet boundary, and hence involve 
nonlocal constraints on the geometry of the world sheet. 

If we follow a line at fixed a in a general planar diagram, we find a sequence of appearances and 
disappearances of Dirichlet boundaries. Thus each such line is associated with a two state system: 
the Dirichlet boundary is either "on" (solid line) or "off" (dotted line). Thus, in addition to the 
target space and ghost worldsheet fields, we also introduce an Ising spin variable s\ = ±1 at each 
site of the worldsheet lattice. A time link joining two + spins is a bit of Dirichlet boundary, and 
one joining two — spins is a bit of bulk. A time link joining opposite sign spins is a spin flip that 
turns a bit of Dirichlet boundary on or off. The sum over all spin configurations then accomplishes 
the sum over all planar diagrams. 

In practice, the specification of the interacting Ising/target space system involves case by case 
technical details, given for <fP theory in Q and for Yang-Mills in ||. Here we only quote the final 
proposal for the </> 3 theory which is the focus of the rest of this article. The amplitude for the sum 
of all planar diagrams evolving an initial state i to a final state / is given by0 
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expj^E N + i - H)(4 + i - 4)0- -pm- pU) + (i - p?)(pUi + pLi)44] J , 

where P- = (1 + s|)/2 and we have put e = me/a. The states \i), \f) are specified by the number 
and position of Dirichlet boundaries that extend to r = — oo,+oo respectively. Here the vertex 
functions are given by 
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§This formula is Eq. (26) of Ref. itself a refinement of the one given originally in Q. It is written with a new, 
more transparent, arrangement of the dependence on the lattice constants. 
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Also the projectors multiplying the vertex functions are 

v\ ee pi pi +i (i - pt'nzM-i n = piPt\\ - pj +1 )ni 1 n| +1; 

where U{ = (1 + s{s{ +1 )/2. 

Without going into great detail we make a few explanatory comments about this formula. First, 
notice that the spin projectors P- keep track of the distinction between solid (P = 1) and dotted 
(P = 0) lines. A bit of Dirichlet boundary is represented by a delta function, identifying the q's at 
successive sites, defined through the limit 

5(Aq) = lim ( ° % e - aAq2 / 2me . (7) 

The prefactors in this formula are taken care of by the ghost integration. Since formally 5(0) = 
Vj_/(2tt) , the volume of transverse space, we should regard e as a temporary infrared cutoff on the 
transverse coordinates, related to the size of the space L = V 1/D by e = 2vra/mL 2 . The different 
arrangement of projectors for ghosts compared to those for q is due to the fact that b = c = on 
solid lines, whereas the q's are only required to be equal on solid lines. The apparently elaborate 
ghost dependence of the vertex functions in (||) provides necessary factors of l/p + that are dictated 
by the field theoretic Feynman rules. The noteworthy feature is that in spite of their complexity, 
they are described locally on the world sheet. The rather complicated combination of projectors 
multiplying these vertex functions limits their occurrence to the appearance or disappearance of 
Dirichlet boundaries. The extra II projectors remove difficulties, due to the ghost insertions in (|5|), 
that occur when two or more vertices are within one or two lattice sites of each other. 

Finally we comment on the limits involved in recovering standard field theoretic perturbation 
theory from Eq. ^. One first expands the formula in a power series in g, T = ^2 n g n T n (e,a,m). 
Then the limit e — > on each T n converts it to a standard Feynman integral, with uv cutoff a and 
p + cutoff m. The dependence on these cutoffs then just parameterizes the standard field theoretic 
divergences. In this article we are attempting to learn about physics at finite coupling by studying 
d|]) in the presence of all these cutoffs. But to compare our results to perturbation theory we must 
keep in mind that the perturbation expansion is carried out before sending e — > 0, implying the 
parametric regime g <C £• 



3 The Mean Field Approximation on a Continuous Worldsheet 
with Simple Cutoff 

In this section, we will develop a somewhat crude, but qualitatively transparent, version of the 
mean field method on the world sheet. A more precise and detailed version of the same method, 
based on the worldsheet lattice, will be presented in the next section. Here we will adopt the 
continuum version of the world sheet, and we will use the notation and conventions of the previous 
section, except that we use the Minkowski metric in both target space (so the factor at each vertex 
is ig, g =coupling const.) and on the world sheet. The worldsheet coordinate a is compactified 
in the interval to p + , and r runs over the interval Tj to r/, which will eventually go to infinity. 
The target space fields q(<r, r), which will be designated as matter fields, live in D dimensions. For 
simplicity, we take the total transverse momentum flowing through the worldsheet to be zero: 

f P daq'(a,r) =q(p+r) -q(0,r) = 0. 
Jo 
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To keep track of the two types of lines, it is convenient to define a scalar field <fi on the worldsheet 
which is 1 on solid lines and on dotted lines. Now consider the expression 

(^) D/2 exp(-W(Vf) (8) 

As (3 tends to oo, this expression tends to a delta function, 5(V), for all 4> 0. Here, we assume a 
Euclidean continuation to have a well defined Gaussian. We will use this to impose the condition 

q = 

on solid lines, where <p = 1. On dotted lines, where <f> = 0, there is no constraint. The matter part 
of the action, with the above constraint, takes the form 

rp+ rr f 

S q = da drC q , 

Jo Jn 

C q = ^>(q) 2 -^(q') 2 . (9) 

For the time being, we will keep the parameter f3 finite, though eventually we will let (3 — > oo. 
Comparing to the lattice expression (|j), we see that <f> here is playing the role of the projectors 
P?P?~ , and = a? /m?e. Also, we have dropped the prefactors in front of the exponential in (JsJ) ; 
these give rise to quantum effects which will not play a role in the approximate treatment of this 
section. However, their effect will be included in the lattice calculations of Sec. 4, where we shall 
see that they are mostly canceled by corresponding factors in the ghost sector. 

In addition to the matter part of the action, there is the ghost action. The function of ghosts is 
to cancel the contribution of the matter fields on the dotted lines, and leave it untouched on solid 
lines. The ghost action is somewhat complicated, and its lattice version will be given in the next 
section. For our purposes, we will not need the explicit form of the ghost action. In fact, we find it 
easier to compute directly the effective actions W q ^{4>) for matter and ghost sectors as a function 
of 4>. The calculation proceeds as follows: The world sheet can be thought of as a union of black 
regions, composed of only the solid lines, and the white regions, composed of only the dotted lines 
(See Fig.l). In the white regions, where = 0, the situation is very simple: The matter and the 
ghost contributions cancel, and the total effective action is zero: 

W q ((f> = 0) + W g (<i> = 0) = 0. (10) 

On the other hand, in the black regions, where <j) = 1, the ghosts do not contribute: 

W g (<t> = 1) = 0. (11) 

Therefore, we need only compute W q (4> = 1) for an arbitrarily shaped black region. This is a 
standard Casimir type calculation. Since the answer is divergent, a suitable cutoff has to be 
introduced, and for the Casimir effect, one has to calculate the cutoff independent finite part. 
Here, we have a much simpler problem: We are interested in the dominant bulk contribution to the 
action, which is going to be the leading cutoff dependent term. It is also going to be proportional 
to the area of the black region and independent of its shape. We will compute this term for a 
rectangular black region of width L in the a direction and length T in the r direction, taking <p to 
be a constant (a and r independent) in Eq.([|). We will also impose periodic boundary conditions 
in both directions and use a cutoff in the mode number. Defining the operator K by 

k = ^dl - dl, 
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Figure 1 : Regions of Worldsheet 



we have 

W q = ~Tr\og{K). (12) 



W q *- A _, ^» g . (13) 



The cutoff dependent part of the action is easily calculated: 

DLT 

Here S, which has the dimensions of a 2 , is the cutoff parameter. It is proportional to m 2 , where m 
is the discrete unit of a on the world sheet lattice. 

We note that the answer depends only on the area LT of the region. Although the calculation 
was done for a rectangular region, it is clear that the leading cutoff dependent term is proportional 
to the area even for a region of an arbitrary shape, and independent of the boundary conditions 
imposed. This dependence on the area, and also the dependence on the combination /3(j>, can be 
established essentially by dimensional and scaling arguments, independent of any detailed calcu- 
lation. The numerical factor 4tt is unimportant and it could be absorbed into the definition of 
the cutoff. For small sized regions perimeter and shape dependent contributions may also become 
important. However, in this section, we will focus on the bulk contribution only, and we will neglect 
corrections of this type. The cutoff 5 hides our ignorance of the texture of the world sheet when the 
continuum limit is taken; it will eventually be absorbed into the definition of the coupling constant. 
Finally, the matter contribution to the action in the black region is obtained by setting = 1 in 
Eq.@: 

DA 

4vr(/3) 1 /2 ( r 

where A is the area of the region. Eqs.(p!o|),([TT|) and ( |l4|) provide us with complete information 
about the contributions of the matter and ghost actions in the two types of regions. 
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Next, we have to tackle the functional integration over (ft. We impose the constraint that (ft = 0, 1 
by means of a Lagrange multiplier 7r(a, r): 



AS'n = / da dr A£ n , 
JO Jn 

AC n = 7r^-2(-l) m 5(r-r m ((7o))^. (15) 

The functional integration over (ft consists of integration over T m {ao) and summation over n. These 
summations and integrations cannot be carried out in closed form; however, the problem can be 
reduced to the solution of a differential equation. Let us set 

exp (ASU)=F i >t(T i ,T f ), 

where AS l ' J is the resulting action after the sums and integrals are carried out. The outcome is 
one of the four possible functions on the right hand side, depending on the boundary conditions on 
(ft at Tj and tj. The indices i and j can take on the values (+) or (— ), corresponding to (ft = 1, at 
r = Tij. These functions are defined through the differential equations 

<9 t F + '+(t;,t) = ig exp(iir(T)) F + ~(Ti,T), 

d T F + ~(Ti,T) = ig exp(-ivr(r)) F + ' + (ri,r), 

d T F~' + (Ti,T) = ig exp(wr(r)) F~~(Ti,r), 

d T F~'~ (tj , t) = ig exp(-ivr(r)) F~' + (Ti,r), (16) 

plus the initial conditions 

F+>+(Ti=T f ) = F-'-(Ti = Tf)=0, 

F + >-( n = T f ) = F-' + { n =T f ) = i. (17) 

In writing these equations, we have assumed that the vertex that converts a dotted line into a solid 
line or vice versa is simply given by the coupling constant g. It was shown in M that there are 
additional contributions to the vertex involving the ghosts b and c; however, in the leading large D 
limit, which will be the approximation scheme adopted in this article, they do not contribute. This 
is because the function of these vertex ghost insertions is to produce the factor of l/p + , which is 
part of the propagator, and it does this by deleting a pair of ghost fields b and c, as explained in 
|l]]. The crucial point is that it is always a pair of fields irrespective of the transverse dimension D, 
since p + always appears to the first power and not to the power D. Therefore, this contribution 
goes like 1 and not D, and it is negligible in the large D limit. We will have more to say about the 
large D limit later. 

It turns out to be convenient, although not essential, to keep track of the (+) and (— ) indices by 
means of two fermionic (anticommuting) variables ei 2{r,a). These fermions, already introduced 
in Q, give the continuum formulation of the Ising system of the next section. e\ is associated with 
the (+) index, or the solid lines, and e2 with the (— ) index, or the dotted lines. Their r dependence 
is given by 

ei(r,cr) = F + ' + (Ti,T)e 2 (Ti,a) + F-> + (T i ,T)e 1 (T i ,a), 

e 2 {r,a) = F + ~{r u T)e 2 {T i: a) + F~'~(n, r)ei{r u a). (18) 
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As a result, the e's satisfy the differential equations 



iei+g exp(?7r(r))e2 = 0, 
ie.2+g exp(— i7r(r))ei = 0, (19) 

leading to the action 

rp+ i-Tf 

Sf = da dr {ie\e\ + ie2&2 + 9 exp(z7r)eie2 + g exp(— m)e2ei) ■ (20) 
For later application, it is convenient to make the change of variables, 

ei = exp ( ~ir\ fa, e% = exp 

e 2 = exp ^2, e 2 = exp (~ir) fa, (21) 

with the corresponding action 

S f = J Q da J dr ^i^i + ^2^2 - ^Trfafa + ^nfafa + 9 tpifa + gfaipi^j ■ (22) 

Adding all the contributions, the full expression for 

T fi = ( Pf \e-^\ Pi ) 

is 



J DqDXDcpDbDc 

(23) 




Tfi 



where C q and Cf are given by Eqs.(|) and ©. 

We are interested in the ground state of the system described by the action given above. We will 
compute the ground state of the combined matter, ghost and the fermionic systems in the presence 
of constant fixed background fields 4> and it, and then minimize the total energy with respect to 
these background fields by solving the classical equations of motion. This is in essence the mean 
field method. A systematic way to do this is to consider the large D limit, where D is the number 
of transverse dimensions. In practice, this number is not particularly large, so the method could 
at best be expected to give qualitative results. However, it is a convenient way of organizing a 
systematic expansion scheme in inverse powers of D. We note that the contributions of the matter 
and ghost fields to the ground state energy are proportional to D, so it is convenient to scale the 
field 7r and the coupling constant g by 

vr -» Dtt, g -» Dg, (24) 

in order to have an effective action proportional to D. In what follows, we will simplify things by 
considering only the leading contribution in D. This does not mean that the non- leading contri- 
butions are unimportant; for example, it is clear that Lorentz invariance can only be understood 
by taking into account the non-leading terms. Also, there is the question of the meaning of the 
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coupling constant g. From Eq.fllq), one sees that this constant always has the dimensions of mass 
(or energy). The dimension of the constant g that appears in the g 3 interaction in field theory 
is mass^ 4-15 )/ 2 , depending on the dimensionality of space-time; only for D = 2 (4 dimensional 
space-time), does it have the dimensions of mass. As a result, we cannot identify the constant 
g that we are using with the field theory coupling constant, except perhaps in four dimensional 
space-time. In other dimensions, we have to treat it as an effective parameter, and to establish 
its connection with the parameters in field theory would again require a knowledge of higher order 
terms in the 1/D expansion On the light-cone worldsheet a has dimensions of momentum and r the 
dimensions of time. On our worldsheet lattice, the lattice constants in the two directions therefore 
have a dimensionful ratio m/a, with dimensions of force. Since this ratio may be kept fixed in the 
limit of a continuous world sheet, the formalism provides a scale which can be used to define a 
dimensionless coupling constant, g = g\J\ /32tt' 2 (a / m) ( 4 --°)/ 4 _ The coupling constant used in ( |l~6| ) 
can therefore be identified with g(m/a)^~ D ^ 4: up to dimensionless factors. 

Another important simplification follows from the nature of the ground state. Since the problem 
is translationally invariant in both the compactified a and the r direction, at least in the limit 
Tf — n — > oo, we expect the ground state to share these symmetries; therefore, we take 4>o, the 
expectation value of <fi, to be independent of a and r. We note that the value of the product (3cj>o 
is especially significant: From Eq.(^), it follows that a finite non-zero value of this product leads 
to the standard string action with a finite slope parameter. The key idea behind our computation 
is self consistency. Starting with a finite (34>o, we compute the zero point energy of the resulting 
string, and add this to the energy of the fermions to get the total energy. Minimizing this energy, 
we arrive at a finite value of (3(po, completing the cycle of self consistency. 

In trying to compute the energy of the combined matter and ghost system for a background 
value (fro of 4>, we encounter a problem. We know the effective action for this system, from which 
the ground state energy is easily deduced, only for <j) = and <j> = 1 (Eqs.©,©, and (0)). On 
the other hand, </>o can take on any value between and 1, so we have to extend the definition of 
the action to an arbitrary 0o between these limits. This can be done as follows: 0o is the classical 
expectation value, or the average value of </>. Consider a specific partitioning of the world sheet 
between the white and black regions, such as represented by figure 1. If we denote the total area of 
the black regions by Ab and the total area of the world sheet by ^4^, and remembering that = 
on the white regions and eft = 1 on the dark regions, the average value of 0, 0o, for this partitioning 
is given by 

00 = ^p. (25) 

On the other hand, from Eqs.(|l^),(|ll|) and (|l4|), the contribution of this partitioning to the com- 
bined matter and ghost action is 

w « +w ° = -imk- < 26 > 

From these equations, it follows that 

W +W - DAw 0° (27) 

which is the equation that expresses the combined action in terms of 4>o- It is now easy to convert 
this into an equation for the corresponding ground state energies, which we label as sj 0) and Ef ] 
respectively. Remembering that the area of the world sheet is given by 

Aw = P + (Tf - n), 
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we have, 

The combined matter and ghost energy therefore has a simple linear dependence on </>q. This is a 
direct consequence of the area law for the matter action (see Eq.([I3|)). Since this simple dependence 
on area is bound to get corrections for small regions, we expect to have some deviation from the 
linear dependence of the energy on <pQ. In fact, a calculation based on the worldsheet lattice, 
presented in the next section, results in a more complicated dependence, but the linear dependence 
can be regarded reasonable approximation. 

Next, we have to introduce some background for the field tt. At first, one might think that, 
by translation invariance, this should again be a constant. However, it turns out that a constant 
background is trivial; it is clear from Eq.(|22|) that only the derivative of tt with respect to r has 
dynamical significance. So we ansatz 

(vr) ^vr (0) (r) =7r r + vri, (29) 

where 7ro,i are constants. Again, we see from Eq.(p2|) that this background, which at first looks 
time dependent, is in fact static, and the system has a well defined energy. To compute this energy, 
it is convenient to quantize the fermion action of Eq.(|22|), in the background given by Eq. (p9[) , and 
construct the Hamiltonian. The result is 

rP + (\ _ l _ 

Hf = da ^-^0^2*02 - -TToV'lV'l - #1^2 - gfalpl 

(o 71 " a n,2 a «-2 - ^7T ali,l a n,l ~ 9 a\±0-n,2 ~ 9 4i,2 a n,l) ■ (30) 

In this equation, a n ^ and a) n i are fermionic operators that satisfy the usual anticommutation 
relations 

where i,j = 1, 2. The fermionic field has the mode expansion 

& = ^(P + y 1/2 a n,i exp (j^^j , 
-1/2 f ( 27rincr\ 



6 = E(p + )" l/a <i«p(-^r-j- ( 31 ) 

n \ P / 

The vacuum is as usual annihilated by the o's. In the state under consideration, each mode number 
n is filled by one of the a^ n j's. This corresponds, according to the value of i, to either having a solid 
or a dotted line. There is no meaning to an unoccupied state. 

The Hamiltonian of Eq.(|30"|) is easily diagonalized. There are two eigenvalues for each mode n: 

4=4 (7r ° 2+V)1/2 - (32) 

To find the energies corresponding to the original fermionic variables e\ t 2, we have to perform the 
transformation of Eq.(21), which results in four possible energies for each mode: 

#i, ± =4^ o 4 (7r ° +V)1/2 - (33) 
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These four different energies correspond to four different boundary conditions which can be imposed 
on the fermionic state at the initial time. However, we shall see that only (— ,+) and (— ,— ) 
combinations satisfy the constraint 

< (f> < 1. 

We will consider both of these possibilities in what follows. 

The next step is the computation of the total fermionic energy, which, as it stands, will be 
divergent. One way to regularize it is to discretize the a interval, with a spacing Act between two 
adjacent points. It is then easy to see that, to go from the constraints of Eq.(|l5|) written in discrete 
form, to an integral over a, we need a factor of Act, which can be generated by a suitable scaling 
of ttq and g: 

ir ^AaiTo, g^Aag. 

Alternatively, one can use a cutoff v on the mode number, the relation between the two cutoffs 
being 

Act — > ttv. 

Using this latter cutoff, the total fermionic energy is given by 

Ef = % £ (-0 - (-o 2 + V) 1/2 ) exp (-^) - P y (-0 + (vro 2 + , (34) 



where, at the end, we have taken the limit of small v 

The total ground state energy is the 
last term in the exponential in Eq.(|23|): 



The total ground state energy is then the sum of Eq°\ Eg°\ Ej and the contribution from the 



E o =DV + (^yr s + -000 - -TTo ± + V) 1/2 j • (35) 

To find the ground state energy, we have to minimize E^ with respect to ttq and Minimizing 
it with respect to ttq gives 

1 dEt 1 , 7T 

- --± — — g + O = O, (36) 



Dp+ Ottq ' 2 2(vr2 + V)V2 

We note that, as ttq ranges between — oo and +oo, the range of 4>o is between and 1, as it should 
be. We plot ^(^o) m Fig- @- Had we chosen the (+,+) or (+,— ) combinations in Eq. (|33"1) , we 
would not have the correct range for <J)q. Next, minimizing the energy with respect to 4>o gives 

1 dEt 1 . , 

" - +^o = 0. (37) 



Dp+ d(p 4vr(/3) 1 /2^ 

We will first consider these equations in the weak coupling limit and compare the results with 
perturbation theory. Accordingly, we will assume g to be much smaller than ttq, and expand Eq,(j35l) 
in lowest order in their ratio: 

= \ ± Q - 16vrV^ 2 ) • (38) 
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-2 -1.5 -1 -0.5 0.5 1 1.5 2 

Figure 2: The mean field <po as a function ttq, as given by Eq. (|3^) with the choice of — sign, for 
various couplings. 



Since in the weak coupling limit we expect the dotted lines (white regions) to dominate, (fio should 
be small compared to one, and the minus sign is the correct choice in the above equation: 

<t) {weak) ^ 16vr 2 c/g/3, (39) 

where we have scaled the coupling constant according to 

9o = g5. 

This corresponds to the combination (— , — ) in Eq.(]33|), and from Eq,(|35|), we see that it has lower 
energy compared to the alternative combination (—,+). It is also clear from Eq.(^) that go must 
satisfy the condition 

4^o (/?) 1/2 « 1 (40) 

for (j)Q to be small, and also, as was assumed, for g to be smaller than Ao- For the sake of com- 
pleteness, we note that the combination (—,+), which has higher energy and is therefore unstable, 
corresponds to (fio ~ 1 and a resulting predominance of solid lines. 

We do not expect perturbative field theory to lead to string formation, so it is important to 
verify this in the weak coupling regime. As we have remarked earlier, strings form only if the 
parameter <f>o(3 is non-zero and finite as f3 — > oo. If we define the weak coupling limit by the 
condition 

5o/3 -» 0, 

then, from Eq.(|39|), we see that there is no string formation. In fact, it follows from Eq.(^[) that, 
as expected, this is the zero slope or the field theory limit. 

The next question is if or when string formation is realized. This happens, at least in the 
mean field approximation, in a regime which we will call the intermediate coupling regime, go still 
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satisfies the condition (39), but now if we require that 



900 ¥= 0, 

which means 

P<fo + o. 

This is the condition for a finite non-zero slope, and quantization of S q then results in the usual 
string spectrum with the level spacing given by ir/((/3(j)o) 1 / 2 p + ). In passing, we note that the 
presence of p + in the denominator is consistent with Lorentz invariance. Remembering that we 
are calculating the spectrum of p~ with p set equal to zero, the product p + p~ has to be a Lorentz 
invariant pure number, as is the case here. Of course, all of our results are conditional on the 
validity of the mean field approximation. 

Finally, we would like to discuss briefly the strong coupling regime, by which we mean the 



range of values of go that do not satisfy the condition given by (40). For example, go could remain 
finite and non-zero as (3 — > oo. This corresponds to the infinite slope limit of the string: The level 
splitting goes to zero, and the whole string spectrum collapses into a continuum. 

The last topic we wish to discuss is the ratio of the density of dotted lines to the density of the 
solid lines. To compute this ratio, which we call r, we need the ground state of Ht. It is given by 



n( 



u„ a{^ + v n 4 j2 ) |0), (41) 



where the constants u n and v n satisfy 

Vn_ 1_ 

u n 2g 
The ratio r for state is then given by 



7^ + V) A/ +to , Kr + M =1- (42) 



\(s\K + datl>Ms)\ Kl 2 ±9 2K 



Since ttq is negative (Eq.(|37D), this ratio is always greater than one, so there are always more dotted 
lines than solid lines. In the weak and intermediate coupling regimes, we can approximately write 

' (44) 



In the weak coupling limit, go/3 — > 0, so r — > oo, even for a finite (3. As expected, this means 
that the worldsheet is taken over by the dotted lines, and there are hardly any solid lines. In 
the intermediate coupling regime, g\j3 is finite and non-zero, so r is finite for finite (5, although as 
(3 — > oo, it is still true that r — * oo. So the dotted lines are still in preponderance, although not 
to the same extent as in the case of weak coupling. It is somewhat surprising that, even in the 
presence of this relative scarcity of solid lines, string formation can take place. 

4 Mean Field Path Integral on the Worldsheet Lattice 

We would now like to apply the mean field method directly to the lattice system (||). It is this 
system which we have shown is equivalent to the summation of planar diagrams O. In particular 
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the worldsheet lattice spacings a, m are directly related to cutoffs in x + and p + that were employed 
to render all Feynman integrals finite. 

In applying mean field theory to Eq. |||, we make the same simplification as in Section 3 by 
dropping the extra ghost insertions at the ends of solid lines (i.e. by setting the exponential factors 
in (||) to unity). As explained in Section 3, mean field theory is expected to be valid when D — > oo 
where D is the number of q fields or twice the number of ghost b, c pairs. In such a limit, this 
simplification should have negligible effect at D = oo, since the dropped ghost insertions involve 
only one of the infinite number of ghost pairs. Once this simplification is made, the extra II 
projectors in (Iff) can be safely removed, and we can then make the replacements 



Because of the projectors, the interaction terms appear only when the ghosts multiplying them are 
decoupled from all other ghosts. This means we can remove the factors 



Am^flT V GOT V 8TO7T 



m /2m£ WT^-V^' (46) 



3 J~ x \ 



from the exponent and then regain the correct answer by inserting the factor Ylij(9') **** m 
the path integral. Here we have also defined the dimensionless couplings g and g' . 
These manipulations lead to the simplified expression 



T fT = S E J DcDbDdl[(2g^) { 



s J =±i " hi 



cxp 



2m ^ 



(qi + i-q|) 2 + (q]-qr 1 ) 2 



PJP: 
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2me 



expj^E [fe " *(4fi - 4)0- - Pi)0- - PL) + (1 - Pl){P! + x + PLi)H4\ |(47) 

Note that the factor yfi multiplying g is a direct consequence of the mismatch of Dirichlet enforcing 
delta functions on each solid line segment: the value of q on the solid line is integrated, whereas 
the ghosts are set to zero everywhere on the solid line. Therefore, the prefactors in the Gaussian 
approximation to the delta functions do not quite cancel. 

To facilitate a mean field treatment we introduce unity in the form 

1 j I k>\[ <)[<,,, /'/) J D^DX exp J i £ Ay - P> ) 1 (48) 

into the functional integrand, and then assume we can treat the A and (j) integrals classically. That 
is, in this section we identify the mean field approximation with a saddle point evaluation of the 
integrals over eft and A. One advantage of this way of implementing the approximation is that it 
sets the stage for a systematic exploration of fluctuation effects that certainly correct, and in some 
cases (as for example with the two dimensional Ising model) invalidate some of the results of the 
mean field approximation. 

The classical treatment of the A integration introduces some freedom in what we take as the 
zeroth order in the approximation scheme. It is this integration that constrains the values of faj to 
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be or 1, and approximating the A integration will relax these constraints to some degree. Thus 
valid rearrangements of the defining path integral (f47[), which exploit the fact that the P's are 
projectors, will lead to different classical (i.e. zeroth order) actions. The fundamental premise of 
the mean field approximation is that slowly varying, and in particular uniform, field configurations 
capture the important physics of the system. If this is really true, the results should be insensitive 
to reasonable rearrangements along these lines. If there is sensitivity, comparing the results of 
different starting points can give some measure of the credibility of the approximation scheme. The 
most naive setup is to blindly substitute (pij for each P? with no further rearrangement. 

How can we decide on the best zeroth order action? Of course, we would like to choose the one 
which gives the closest possible agreement with the actual answer. For coupling the mean field to 
the matter and ghosts, we can get some guidance from a simple exact evaluation of Eq. |] at zero 
coupling. Then all solid lines are eternal, and we can easily calculate the exact energy for n equally 
spaced solid lines that all extend from early to late times. For n/M fixed this should correspond 
to the energy at uniform mean field (f) = n/M. By construction, the formula gives for the integral 
over all variables on dotted lines 

N / n \nD/2 ( _ n 13 _ n 3\2 n ) 

n(^) ^{-^^J^-^H-^?}, («) 

where these q's are those on the solid lines. Integrating over all of them gives a lattice string path 
integral of the type exactly evaluated in § and briefly discussed in the appendix. The answer for 
the bulk energy is 

= M </>Re{iLi 2 (-iv^e)} (50) 

„ M —^l\ (51) 

where in the first line we have taken M large, in the second line identified n = M(j>, and in the 
last line taken e small. There are of course other configurations of eternal solid lines that could 
simulate a uniform mean field, but this one has the lowest energy. 

As far as mean field calculations go, given their relative crudeness, taking the exact zero coupling 
result (|5^) as the matter and ghost contribution to the energy, or even taking the linear interpolation 
used in Section 3, is probably as reliable as other choices. However, if we want to use the saddle 
point formulation to systematically go beyond this approximation, it is useful to have a zeroth 
order action that leads to a qualitatively similar answer. Also, we would like to identify the mean 
field in some way with an effective string tension, and for this we need to know at least how the 
target space fields couple to the mean field in the path integral action. Finally, once interactions 
come into play, the lowest zero coupling energy given by (]50| ) need no longer be the appropriate 
energy to assign to the ghost and matter system, and other possibilities should be kept in mind. 

Even so, the result (|50| ) is strong support for immediately disposing of the naive choice for 
zeroth order action, which leads to the bulk ghost-matter energy density 

£ » r = ^ { !Ll2 (f)}- D >„ + , (52) 
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Figure 3: Comparison of the ghost+matter energy for the naive zeroth order action (dashed lines) 
to the zero coupling result of Eq. ^ (solid line) at e = 0.01. The higher dashed curve has the ghost 
mass parameter fi set to zero. 



where [i = e(l — <fi)(2 — (f>)/2. In Fig ||, we plot on the same graph this result, the same expression 
with \i = 0, and the exact zero coupling one. We see that the curves for (jaC]) and (|5^) disagree 
in almost every respect qualitatively and quantitatively. Clearly, we should find a better starting 
point! 

One problem with the naive starting point is that the non-gradient ghost terms do not have a 
transparent formal continuum limit. This bodes ill for a slowly varying mean field approximation. 
We can improve this situation with two rearrangements. First consider the change of ghost variables 
b{, c{ — > bi — -P/&| _1 ,cj — -P/c- -1 , which has unit Jacobian. Then we see that the only change in 



Eq. 47 is the replacement of the last term in the second line, 



The shift has no effect on the other ghost terms because they contain at least one factor of 1 — P? . 
Now the formal continuum limit of ( |53| ) is transparent. As stressed in Section 3, the function of 
the ghosts is to cancel the effects of the matter fields on dotted lines (where fyj = 0). This change 
makes the cancelation more effective for small but nonzero mean field as well. We shall refer to the 
ghosts treated this way as dynamical ghosts, and shall restrict the calculations in this section to 
that case. Secondly, we can rearrange the last term of the last exponent in ([!?]) which as it stands 
also has an obscure formal continuum limit. Note the identity 

(i - pixpUi +p L) = \ [ p U +pu-zpi+(pi- pUif +( p i- ^ti) 2 ] • ( 54 ) 

Although it takes more space to write, the right side is formally of order m 2 , so the formal continuum 
limit is apparent. In addition this rearrangement causes this "ghost mass term" to vanish for 
uniform mean fields. In all that follows, we shall adopt both of these rearrangements. 
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If we make only the rearrangements mentioned thus far, substitution of (f>ij = (f> leads to a 
ghost and matter path integral action in which acj) 2 /me multiplies the time difference terms of 
both matter and ghosts, and a/m (a(l — </>) 2 /m) multiply the space difference terms of the matter 
(ghosts) . The upshot of doing the path integrals with this action is the bulk energy 

Here we see the basic delicacy with the saddle point approach to the mean field approximation. 
The result agrees with the exact answer at (f> = 1, but the 4> — > behavior depends on a delicate 
cancelation between the matter and ghost contributions. In particular, unlike the exact calculation, 
the limits e — > and (p — > do not commute. 

Another reasonable choice for zeroth order action emerges from rearranging the products of 
projectors with the identity 

Pi Pi - 2 2 2 2 ' 1 j 
and we could similarly write 

pi I pi tpi _ pi\2 

(1 " - P!) = 1 - % 1 - 1 1+1 2 tj • (57) 

If we do both these things we arrive, for uniform mean fields, at the bulk energy 

We compare the different possibilities in Fig |] for e = 1.0 and in Fig. [5] for e = 0.01. 
We see that in Fig. || the curve for a compromise form £ 3 

lies virtually on top of the exact answer. This is because the squaring of in the second term 
makes both curves approach zero like (p 3 / 2 and the squaring of (1 — (f> 2 ) removes the square root 
branch point at eft = 1. The zeroth order action that leads to this interpolation is given by the 
substitution (1 — -P/)(l — -P/+i) —> (1 — ^i 2 )(l — ^i+i) m front of the spatial difference term for 
ghosts, and the substitution (^fj) in front of the time difference terms for both matter and ghosts. 
It thus involves a valid but non-minimal interpretation of the projectors appearing here. Because 
it does well with the zero coupling case, we shall take this zeroth order action in the mean field 
calculations to follow. 

Recall that the saddle-point approach to the mean field approximation consists of replacing 
integration over Ajj with the classical equations of motion for these variables. This procedure 
factorizes the remaining quantum averages into the Ising sum times the Gaussian integrals over 
q, b, c. To study the ground state we look for a solution with constant (ftij = (ft and Xij = A. The 
formal continuum limit of the the zeroth order matter and ghost action described in the previous 
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Figure 4: Comparison of ghost+matter energies for e = 1.0. The solid line is the exact zero coupling 
energy (j50j), the highest dashed line is £ gg , the dashed line with infinite slope at (j> = 1 is £^ g , and 
the third dashed line is the compromise . 
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paragraph is for uniform^] mean field 4> 

S q , 9 = J drda (iq* + ^ - (1 - <P 2 )H'c' - ^ (60) 

from which we see that the effective string tension as a function of the mean field is T e //(</>) = 
m-y/e/ay^. In the following we also use /? = a 2 /m 2 e, in terms of which the effective tension reads 
T e ff(4>) = (yT?^) -1 . We are of course interested in the limit f3 — ► do, which, if taken at fixed <f>, 
implies T e ff — * 0. However, the perturbative field theory is characterized by a worldsheet which in 
the bulk has no q 2 term which means T e ff ~^ °°- The perturbative situation will be recovered if 
we find that, as a function of (3, (ft tends to faster than 1//3. 

As a final preliminary to the mean field calculation, we note that the Ising sum can also be 
done exactly for constant mean field A. The sums for each i factorize, so we have 

£ ( ~ (i-^rV2 e -^^= i £ fflt-M-'We^uA . (61) 
s{=±l \si=±l J 

where we put g 1 = 2g^fe The sums over spins just amount to matrix multiplication of the transfer 
matrix 

M e / ) <«> 

by itself iV times. The ground state energy is thus —M/a times the In of the largest eigenvalue of 
this matrix. The eigenvalues are easily found 



1 + e~ iX ± J(l - e~ iX ) 2 + 4g> 2 e~ iX 
t± = , (63) 

with the corresponding unnormalized eigenvectors given by 

We shall find that the classical equations for <p, A imply that iX is real, in which case t+ is the 
largest eigenvalue. We therefore find that the ground state energy of the Ising system is 

aE, = —M 1, 1 + ^(1-^ + 4^ ((j5) 

where we have put — iX = n in anticipation that it is real. It is also of interest to give the ground 
state expectation of (1 ± s)/2 given by 

v it t n( 1 °\ T n v -i 
= WT^Tv = 17^~ (66) 

v' T T N (° °)t n v 
(G\——\G) = hm mriw 1 = 7 r 67 



^The formal continuum limit for slowly varying non-uniform fields has only one new term: 

dTda ( h' 2 + <5^ 2 - (1 - * 2)2b ' c ' ~ 5^ ~ \w + 2 *' 2)bc 
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Figure 6: The mean field </>(/«) of Eq. ( |70| ) for various couplings c = 16g 2 e. Note the similarity to 
Fig. 0. 



These two quantities give the mean number of solid lines and dotted lines respectively. 

• (with i 
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Putting all of the contributions to the ground state energy of the system together (with 



chosen for the matter and ghost contribution), we have 



M 7T 1 \iV4> ) J 7T 1 y iy/(f> l\ 

, , 1 + e K + JU - e K )' 2 + ig r2 e K 
+K(b — In . 

The last step is to minimize with respect to k,,4> and take (5 — > oo. We first minimize with 
respect to k obtaining 

-l + e K + 2g' 2 + y7T^>) 2 + 4^e* 

W v /(l-e K ) 2 + 4^e K 1 + e K + y/(l - e K )' 2 + 4g r2 e K 1 ' 

e K - 1- 1 sinh(K/2) , N 

= = o + — ; = (70) 



2y , sinh 2 (fv/2) + 

which confirms that = (-P/). For g' 2 <C sinh 2 (K/2), — ► or 1 according to whether k < or 
> 0. Thus </>(/«) is a smoothed approximation to 9(k), with the step getting sharper as g' — > 0. In 
Fig. |6] we plot as a function of k for various values of c = 4g' 2 . 

To give an energetic interpretation of the solutions to these equations, we first solve (70) for n 
calling the solution kq{4>): 

0W \2y/<Kl -<(>)} 1 VW " J) J 
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Figure 7: The function — determined from the minimization of the energy with respect to 
for various values of e. We have set D — 2 for this plot. 



Then we substitute Ko((p) for k in the energy to get £{<j), kq (</>)) as a function of <j>. In fact, the 
dependent part of the energy becomes simply 



K(j) — In 



1 + e K + V(l - e K ) 2 + 4g' 2 e K 



<pKo(4>) - In 



- 0) + ^e(l - W + 5^(1 - 20) , 
Minimization with respect to eft gives a formula for k 



<<t>) 



7T<6 



arctan 



1 + 30 2 /^(l-0 2 
=- arctan 



Clearly the stationary points of the energy as a function of <j> are given by 

^-£{<t>,Ko(4>)) = Kq{4>) - k{4>) = 0. 



(72) 



(73) 



(74) 



In Fig. we plot — n for several values of e. (in this and all the plots that follow we set D = 2, 
corresponding to 4 dimensional space-time.) A noteworthy feature here is that —k stays positive 
over the whole range of <f>. As a consequence, the solutions of k = kq will all be in the interval 
< < 1/2. In Fig. |8] we plot £(<p, Ko(</>)) for b = 0.01 and c = 0.1. It shows a single stationary 
point which is a minimum. 

We can analyze the solutions analytically in the limit e — > 0. Solutions with <j> ^ 0, 1 in the 
limit are obtained by expanding k and k$: 



3DJTc 



TT 



K 



2gVe(l - 20) 



(75) 
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Figure 8: The energy curve £(</>, kq(</>)) f° r e = 0.01 and c = 16g 2 e = 0.1 is the lower one. The 
upper curve uses the linear interpolation for the matter+ghost energy as in Section 3. We see that 
the differences between the two treatments are quite mild. We have set D = 2 for this plot. 



The ^fi cancels, and the equation for the minimum becomes 

ZDy/4> 2g(l - 2<p) 



7T " 4>) 



0. (76) 



The left side is plotted for various g values in Fig. ^. As already mentioned, the solutions all lie in 
the range < (j) < 1/2. Notice that Eq. ( |7"6| ) simplifies to a cubic equation for <p with coefficients 
linear in g 2 . Thus the solution 4>{g) as an analytic function of g will have branch points in the finite 
g 2 complex plane, and so its Taylor expansion about g = will have a finite radius of convergence 
as expected for the sum of planar diagrams [||. By virtue of the identification of the string tension 
T e ff((p) = l/y/JJcp = ra^fej{a\f4>\ we see that all solutions for g finite as e — > correspond to zero 
string tension. In other words the system falls apart into a tensionless soup. 

Next we analyze the case for eft near zero. We define a rescaled mean field by r/ = 0/e and 
take e — > at fixed r/. Note that this is the combination that enters the effective string tension 
T e ff = l/y^0 = m/\/ a 2 r], so a solution with finite r\ in the limit e — ► would indicate a finite 
non-zero effective string tension suggesting string formation. Then we find for the energy in this 
limit 

£(ne,K) e(e r? arctan 2r7sinh _1 -^ 2t? f | . (77) 

V 71 Vv Vv g + Vv + g 2 / 

Applying the same limit to k, we find 

De ( 1 if /2 \ , x 

n{eri) « -— I 4t? arctan — - ^— \ , (78) 
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Figure 9: The left side of Eq. fl7q), which determines the mean field at finite g for various couplings. 
Note that all solutions are in the range < <f> < 1/2, with strong coupling corresponding to the 
upper end of this range. We have set D = 2 in this plot. 



and the slope of the energy is given by 

d£(cf), K {(f>)) 



+ K Q ((f>). 



(79) 



The stationary points are of course the solutions of re = kq. Since k = O(e), we can only get a 
solution at fixed r] if we also have g = 0(e). Putting g = Ge, the stationarity condition reads 

(80) 



D A 1 
— 4r/ arctan 

I" V VV 



V 

1+77 



2G 

Vv 



o. 



It is clear from the fact that for any G > 0, the left side goes to — oo (+oo) when r] goes to (oo), 
that there is always a solution to this equation. Note incidentally that had we used ( |50| ) for the 
ghost+matter energy the first term on the left would have been replaced by 3D^/rj/7r, that is, by its 
asymptotic form at large r], and the same qualitative conclusion applies. But there are quantitative 
differences in the results. For comparison of the energy derived from ( |50|) in the fixed r\ regime 
with g = O(e) to ( f77| ) in the same regime we plot both energy curves for G = 1 in Fig. |l^. 

Finally we note that for G <C 1 or g <C t, there is a solution for rj <S 1 given by r] w (G/D) 2 / 3 , 
which corresponds to effective tension T e ff ~ m(D/G) 1 ^ 3 /a. This limit corresponds to the regime 
of perturbation theory where the worldsheet fields in the bulk of the world sheet are constrained 
in terms of the fields on the boundaries. 



5 Conclusions 

In this article, we have presented a scheme for approximately summing Feynman graphs of arbitrary 
order. The scheme is an adaptation of the mean field method, widely used in many body physics 
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Figure 10: Energy at fixed r\ for G = g/e = 1. The solid (lower) curve is from ([77]) and the dashed 
(higher) curve uses ( |50| ) for the matter -fghost energy. For this plot we set D = 2. 



and field theory, to the worldsheet representation of perturbation theory that we have developed 
earlier. For simplicity, we chose cf) 3 as a toy model for study, although we hope to apply the 
methods developed here to more realistic models, such as non-abelian gauge theories. The goal was 
to investigate under what circumstances string formation is possible in field theory. The worldsheet 
reformulation of field theory seems ideally suited to such an investigation; in particular, the field (p, 
which naturally appears in this formulation and which roughly represents the density of propagating 
particles, is closely related to the string tension. A non-zero expectation value for this field signals 
formation of a string with finite tension. 

The problem is to find a reliable way to compute this expectation value. In this article, we have 
used the mean field method as a possible approach to this problem. In this approach, the fields 
on the worldsheet are assigned various background expectation values, which are then calculated 
self consistently by minimizing the total energy of the system. In this fashion, one can find out 
whether the energetics can support a non-zero string tension. Unfortunately, there are many 
technical difficulties in carrying out this program. The biggest problem is connected with the 
delicate cancelation between the matter and ghost sectors, where two big numbers can cancel each 
other to give zero. Because of this, the results of any approximate treatment can be somewhat 
unstable. We try to bypass some of these problems by adopting a simple minded approach in the 
continuum version of the worldsheet in section 3. The calculation supports string formation in 
the intermediate and strong coupling regimes. In section 4, the mean field method is applied to 
the worldsheet lattice, where various details can be treated more exactly than in the continuum 
approach. There the mean field approximation was identified with a certain saddle-point evaluation 
of the worldsheet path integral. Despite its more solid starting point, this approach also suffers from 
possible ambiguities in the choice of zeroth order actions for the approximate evaluation. These 
ambiguities are partially resolved by choosing a zeroth order action that has a transparent formal 
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continuum limit and agrees well with an exact calculation of the contribution of equally spaced solid 
lines in the zero coupling limit. The outcome of the lattice calculation, which also results in string 
formation in the strong coupling regime, is in qualitative agreement with the results of section 
3, although there are some quantitative differences. The most notable of these is the 3/2 power 
dependence of the matter energy on the field 4> (Eq. (|5T|)), as opposed to the linear dependence given 
by (28). We attribute this difference to the different ways of handling the matter- ghost cancelation 
mentioned earlier, and we are moderately encouraged that there is qualitative agreement. We also 
recall the point made in Section 3 that neither the continuum nor the lattice methods adequately 
treat the so-called small black regions. These correspond to fluctuations of the field <f> over regions 
of the size of a small number of lattice sites. The mean field approximation, which at least in the 
leading order, treats </> as a constant over the whole worldsheet, neglects these fluctuations. There 
is indication, for example from the calculation of density of solid lines at the end of section 3, 
that these small black regions may play an important role in string formation. The identification 
made in Section 4 of the approximation with a saddle-point evaluation provides a useful setting to 
systematically explore the ramifications of these fluctuation effects. 

Because of its various shortcomings mentioned above, we consider the dynamical calculation 
presented in this article as a promising initial attempt, not the final word. We hope that the 
line of approach developed here will stimulate further progress in this important problem of string 
formation in field theory. 
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A Determinants 



All of the Gaussian integrals we need in this paper are of the generic form 

1 = /^expj-^^-^V + ^+i-^) 2 ]! (81) 

These integrals can be done exactly using the methods of For simplicity we impose periodic 
boundary conditions in a and Dirichlet conditions in r, which describes a closed string moving in 
a background of Ising spins. 

The eigenvalues of the bilinear forms in the exponent are well-known, so the integral which is 
proportional to one over the square root of the determinant of the bilinear form can be written as 
a double product, so we find for two bosonic variables: 

/ \ MN N M-l , t-> t \ — l 

* = ( T ) n n (* + 4 sin 2 (82) 

W n =i io ^ 2(iV + l) AM) y ' 

The analogous ghost integrals will be the reciprocal of this expression (with of course different 
values for A, B. We list some infinite products from Q: 

n X / ,/ nir \ \ sinfiV + 1)k 9 k , . 

Usin 2 — — A-z\ = — — , z = 4sin 2 - 83) 
V \2(N + 1)J J s inK 2 V J 
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W L • 2 J* \ sin 2 MX A . 2 x , 0/) . 

4 sin 2 — -z = 5—, z = 4sm 2 A (84) 

11 \ MI sin 2 A 



M-l / / \ AT / N 

n (*- | _ m, n (*- jjj^y - to ^ 

Z=l n=l v ' 

the last line being just the case z = of the first two. Using these products, we then find 

* _ J-fir'nf^" 6 )", rf | = ^ 8in ^ (86) 

The ground state energy associated with this system of two bosonic variables can be read off by 
identifying the coefficient of — N in lnA4 as N — > oo: 

where we note that sinh _1 (z) = ln(z + vT-j- z 2 ). So far everything is exact. Now let's consider the 
continuum limit, M — > oo. Note that since a,m — > together in the continuum limit, it is valid to 
keep a/m fixed. We also keep A,B fixed. The bulk contributions to this limit are straightforward: 



° ^ dx sinh 1 f —== sin xir ) — In — = —Re < iLi2 [ ) I — In — 



M Jo \y/A ) A n \ \iVAjj A 

We note the appearance of the dilogarithm or Spence function Li 2 (x) = dilog(l — x) = J2kLi x n /n 2 . 



B Slowly Varying Mean Fields 

The identification of the mean field approximation with a saddle point approximation to the path 
integral, enables an easy extension to slowly varying fields. Indeed we have already seen in the 



footnote to Eq. |60j, the relatively simple extension of the effective action for ghost+ matter fields 
to slowly varying fields. Here we sketch the corresponding extension for the spin system sum. This 
will give dynamical terms to the mean field. 



We first write out the saddle point equation for general (non-uniform) (p, A. 



=±i 

<>kl 



(89) 



This equation implicitly determines Xij{(/>) as a function of all the (fiij. The equation is intractable 
for general <f>ij, but we have seen in the main text how to explicitly solve the equation for uniform 
(ft, obtaining Ay = i«o (</>)• We can get a perturbative solution for slowly varying fcj = 0o + S<f>ij- 
Actually, rather than A itself, we are more interested in expanding the effective action 

W.faXtf)) = »^ /(A) (90) 

ij 

e -ifW = (g'f- s ^ 1)/2 e- l ^ X '^ (91) 

sl=±l 



2(3 



about <pQ. To expand to second order in 8(f), we only need to compute first and second derivatives 
of W s , which by virtue of the fact that W s is a Legendre transform, are simply: 



^ = \A<h) (92) 



d 2 W s d\ l3 



d4>mnd4>ij 

Then the expansion to quadratic order reads 

1 <9A 

W s {<t>, \{<j>)) = W s (<j> , X(<P )) + ^«»«o(^o) + o E ' ' ' ' 



2/ , ~ Y"fc / " y inn 



mn 

mn 



(93) 



+ O(50 3 ). (94) 



The matrix dX/dcf> is the inverse of the matrix d(f>/d\, which can be related to spin correlators by 

-i[{Pip-)-{Pi){p-)\ (95) 



d(f>ij 



dX 



inn 



01) = S * —i ^ 7-. (96) 

i 

Since these spin correlators are all in the presence of a uniform background field A, they may be 
explicitly evaluated in terms of the eigenvalues of the spin transfer matrix introduced in Section 4. 

<«> = w>«> + s lm 4Ba+ £, (iW2)) <m 

so that 

8X mn Wm \t + ) 4(5' 2 + sinh 2 ( K0 /2)) ' 1 j 

The matrix Mj n = p\i~ n \ can be inverted by defining Wj = J2n P^~ n ^v n and proving the recursion 
relation 

Wj+i + Wj-i = (p + l/p)wj + (p - l/p)vj , (99) 

from which it follows that 

(M~ 1 )j n = p2 > _Y ( < ^J+ 1 »" ~*~ ~ 2^jn) ~ ~ ^ #jn 



(l-j' 2 )(W + ^-i.-2g , ^ yg /2 + S mh 2 ( K /2) ^ 
4cosh(K/2) v /5' 2 + sinh 2 ( K /2) ^ cosh( K /2) 

Plugging these results into the quadratic approximation to the spin system effective action then 
gives 

iw s (4>, x(4>)) « iW s {<i> G , inoifo)) - ]T WiMfo) - \ E " ^,i-i) 2 + /^ 2 ^] (ioi) 

ij ij 

(l-^ 2 )V / ^ 2 +sinh 2 (K/2) 2 _ (^+sinh 2 (K/2))3/ 2 

5' 2 cosh( K /2) ' ^ ~ <?' 2 cos1i(k/2) ' 1 j 
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This result gives the spin system effective action for <fiij = (f>o + 54>ij to quadratic order in 5(f). 
Formally <pQ can be chosen to be anything, but to minimize corrections to the slowly varying 
mean field approximation, it should be chosen so that Ylij S&ij = 0- I n the uniform mean field 
approximation studied in the main text, this means simply = (fi = cft . 



It is instructive to combine the work of this appendix with that in the footnote to Eq. |60j to 
give the formal continuum limit of the total effective action to order 5(f) 2 : 

S = -iW.fo.iKofo)) + / d ^<? (U 2 + T^q 2 - (1 - " ^rbc 

J \2 2m z e m z e 

-h&" + 20' 2 )6c + ^-(f) 2 + ^(0 - </>o) 2 ) + 0((f) - (f) f. (103) 

This equation expresses the mean field dynamics for our system in a form tantalizingly close to the 
outcome of the Maldacena conjecture Q. The mean field has emerged as a dynamical Liouville-like 
worldsheet field that enters the worldsheet action in a way strongly reminiscent of the way the radial 
AdS coordinate enters the worldsheet action for the "string" description of M = 4 super- Yang- Mills 
theory. 
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